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We construct a quantum algorithm tiiat creates tiie Laugiilin state for an arbitrary number of particles n in the 
case of filling fraction one. This quantum circuit is efficient since it only uses n(n — l)/2 local qudit gates and 
its depth scales as 2n — 3. We further prove the optimality of the circuit using permutation theory arguments 
and we compute exactly how entanglement develops along the action of each gate. Finally, we discuss its 
experimental feasibility decomposing the qudits and the gates in terms of qubits and two qubit-gates as well as 
the generalization to arbitrary filling fraction. 
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One of the most important goals in the field of quantum 
computation is to achieve a faithful and efficient simulation of 
relevant quantum systems. Feynman first suggested the 
possibility of emulating a quantum system by means of other 
specially designed and controlled quantum systems. Nowa- 
days, ultra-cold atoms in optical lattices are producing the 
first quantum simulators IJ] and many theoretical proposals 
have been already presented in order to approximately simu- 
late strongly correlated systems (0, H 0, lE H] . 

A more ambitious goal consists in finding the exact quan- 
tum circuit that underlies the physics of a given quantum sys- 
tem. Rather than searching for an analogical simulation, such 
an exact quantum circuit would fully reproduce the properties 
of the system under investigation without any approximation 
and for any experimental value of the parameters in the the- 
ory. It is particularly interesting to device new quantum algo- 
rithms for strongly correlated quantum systems of few parti- 
cles. These could become the first non-trivial uses of a small 
size quantum computer. It is worth recalling that strong cor- 
relations are tantamount to a large amount of entanglement in 
the quantum state and this, in turn, implies that the system is 
hard to be simulated numerically ifioll . 

An exact quantum circuit would start from a simple product 
unentangled state and would create faithfully any desired state 
on demand. That is, such a quantum algorithm would diago- 
nalize the dynamics of the target system. At present, very few 
cases are under control ^UM- InRef. Hi, the underlying 
quantum circuit that reproduces the physics of the thoroughly 
studied XY Hamiltonian was obtained. The philosophy in- 
spiring that circuit was to follow the steps of the analytical 
solution of that integrable model. It is not obvious how to de- 
sign a quantum simulator for non-integrable systems. Here, 
we shall present a quantum circuit that allows the controlled 
construction of a particular case of the Laughlin wave func- 
tion. Thus, the quantum algorithm we are putting forward 
will not produce the complete dynamics of a Hamiltonian but 
rather a specific state. That is, our quantum circuit will trans- 
form a trivial product state into a Laughlin state with filling 
fraction one by means of a finite amount of local two-body 
quantum gates. 

Let us start by recalling the Laughlin ifisll wave function for 
filling fraction ly = 1 /m, which corresponds to 
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where zj = xj + iyj,j = l,...,n stands for the position 
of the j-th particle. This state was postulated by Laughlin as 
the ground state of the fractional quantum Hall effect (FQHE) 
ifTill . From the quantum information point of view, the Laugh- 
lin state exhibits a considerable von Neumann entropy be- 
tween for any of its possible partitions ifTsll . It is, thus, classi- 
cally hard to simulate such a wave function making it an ideal 
problem for a quantum computer. 

We shall construct the Laughlin state using a quantum sys- 
tem that consists of a chain of n qudits (d-dimensional Hilbert 
spaces). In our case, that is to = 1, the dimension of the qu- 
dits is needed to be d = n. Let us proceed to construct the 
quantum circuit by first considering the case of = 2 parti- 
cles, then n = 3 and, finally, the general case. 

The Laughlin state can be written in terms of the single par- 
ticle angular momentum eigenstates, also called Fock-Darwin 
states ifi{z) = {z\l) = exp(-|z|2/2)/^/7^I!. Then, the 
n ~ 2 Laughlin state reads 



1 



V2 



{ipi{zi)ipo{z2) - ipo{zi)(pi{z2)) ■ (2) 



Let us note that Z2) = (zi, Z2|^'2^^) 

projection of the Laughlin state 
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is simply the 
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in coordinates representation, where particle label is retained 
in the order of qubits and the angular momentum or 1 is an 
element of the angular momentum basis. 

It is trivial to find a quantum circuit that transforms a prod- 
uct state into the above n = 2 Laughlin state. Let us first 



prepare an initial state as j'l'g 
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|0, 1) and perform on it 



the simple two-qubit gate U[2] as shown in Fig. [T] The exact 
form of the unitary operator U[2] in the angular momentum 

basis{|0,0),|0,l),|l,0),|l,l)}is 



U[2] = 
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Let us now move to the more complicate case of the Laugh- 
lin state with three particles, n = 3. We now need to consider 
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FIG. 1: Scheme for the quantum circuit that generates the n 
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10, 1). Notice 



Laughhn state by acting on the product state |*I'q 
that the qubits of the input state are labeled from bottom to top. 
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acting on a product state Iv&o'''') = jO, 1, 2). 



3. Following similar steps 

(3)\ _ 



a system of three qutrits d = 

as we did for n = 2, we take as initial state I^'q'''') = |0, 1, 2) 
that is, each qutrit is prepared in a different basis element, rep- 
resenting different angular momenta. The aim of the quantum 
circuit is to antisymmetrize this initial state, since the Laugh- 
lin wave function for m = 1 is simply the Slater determinant 
of the single particle wave functions 
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To do this, we define the two-qutrit unitary operators Wjj {p) 
as 



W^J{p)\ij) = ^/p\ii) - ^/l-p\ji) 



(6) 



for i < j, < p < 1, and W^j\kl) = \kl) if {k,l) ^ 
We realize that for the case of qubits and p = 1/2 we recover 
the gate U[2] of Eq. (|4|i. Let us note that the unitary oper- 
ator Wij is a linear combination of the identity (p ~ 1) and 
the simple transposition {p = 0) operators, where a simple 
transposition is defined as the transposition between two con- 
tiguous elements. The architecture of the quantum circuit that 
produces the 71 = 3 Laughlin state in Eq. (|5]l by means of the 
local gates from Eq. (|6]l is presented in Fig. |2] 

So far, we have seen the quantum circuits that produce the 
Laughlin state for n ~ 2 and 3. From these cases, a general 
scheme emerges that will produce the correct quantum circuit 
for an arbitrary number of particles. Let us proceed by induc- 
tion. We will assume that we already know the quantum cir- 
cuit, U [n] , that produces the Laughlin state for n qudits when 

acting on I^'q"'') = |0, 1, 2, . . . , ?i — 1). We now need to com- 
plete the circuit to achieve the n + 1 Laughlin state from the 



product state | 



= |o,i,. 



l,n). 



The Laughlin state for n qudits has the form 
^^sign(7')|ai,...,a 
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FIG. 3: Quantum circuit that produces the Laughlin state for 
an arbitrary number of wires (n + 1) acting on a product state 
, n). Note its recursive structure, I7[n + 1] — 
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FIG. 2: Quantum circuit that produces the n = 3 Laughlin state ''1 ''i 



yI"+ ^1 17 ^ where V gates are defined in Eq. l[9j . 



where the sum runs over all the possible permutations of the 
set {0, 1, . . . , n — 1} and, given a permutation, Uk represents 
its fc-th element. The relative sign between two permuta- 
tions corresponds to the parity of the number of transpositions 
needed to transform one into the other If we add another qu- 
dit, the system is in a product state j^l'^"'') jn). According to 
Eq. (|7]i, we want to generate a superposition of {n + 1)1 per- 
mutations corresponding to the n + l states of the Laughlin 
wave function, from the superposition of ri! permutations that 
we already have in the n-qudit case. 

Let us note that, if we have the set of ??! permutations 
of n elements, we can generate the set of permutations of 
71 + 1 elements by performing successive simple transpo- 
sitions between the new element and its preceding neigh- 
bour in the sequence ({ai, . . . , a„, 71}, {ai, . . . ,n, an}, ■ ■ ., 
{n,ai, . . . ,an})- This idea suggests a circuit as the one 
shown in Fig. [3] In this scheme, the gate V^l"^^' should pro- 
duce a superposition of all the permutations |ai, . . . , a„, ri) 
with \ai, . . . ,n, a„). The gate should do the same task 

in the next site, that is, should produce a superposition be- 
tween |ai, . . . , a„_i, 71, a„) and |ai, . . . , 77, a„_i, an)- This 

scheme works successively till f/"^^'. 

This general structure implies that has to be decom- 

posed in terms of the VK-gates presented previously, 

- Won{p)Win{p) . . . Wn-mip) , 
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where p is a common weight due to the fact that the 0, . . . , rt — 
1 states in the Laughlin wave function are indistinguishable. 
Let us note that all the operators Wm in the previous expres- 
sion commute among themselves and, therefore, the order in 
which they are applied is irrelevant. 

In order to determine the weight p in the V^l"^^' gate, we 
realize that if all the transpositions only involve the state 77, the 
states |ai, . . . , a„, 77) will not be affected by the rest of gates, 
and they should already have the correct normalization factor 



after applying V,i . This implies 
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We proceed in a similar way to determine the rest of gates, 
and we obtain 
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FIG. 4: Quantum circuit tliat produces the Laughlin state of 5 parti- 
cles acting on a product state |01234). 



where, given a k, all the W gates have the same weight -j^^ 
and k = 1 ... n — 1. Notice that the Wij gates act on states 
with i < j along the whole circuit and, therefore, they always 
generate the negative combination — y/1 — p\ji)- This 

is the reason why each term of the final state has the appropri- 
ate sign, since a minus sign is carried in each transposition. 

The above discussion produces our main result, that is, the 
circuit shown in Fig. [3]that uses the definition of its gates in 
Eq. (|9]l. Such a quantum circuit will generate the m = 1 
Laughlin wave function for an arbitrary number of qudits. In 
particular, the quantum circuit corresponding to 5 qudits is 
presented in Fig. |4] 

The recursive structure of the circuit (Fig. [3]l makes easy to 
calculate how the number of gates N{n) and the depth D{n) 
of the circuit scale with the total number of particles n. An 
elementary counting gives the result 



N{n) = {n- l) + N{n- 1) = 
D{n) ^ D{n - 1) + 2 = 2?! - 3 



n{n — 1) 



(10) 



with N{2) = 1 and D{2) = 1. The quantum circuit that de- 
livers the Laughlin state is, thus, efficient since the number of 
gates scales polynomially. This is a non-trivial result since, in 
general, an arbitrary unitary transformation requires an expo- 
nential number of gates to be performed. 

Let us now discuss the optimality of our quantum circuit. 
As we mentioned previously, a simple transposition, s;, is de- 
fined as the transposition between two contiguous elements, 
i and i + 1. Any permutation can be decomposed in terms 
of a series of simple transpositions and its minimal decompo- 
sition is called the canonical reduced decomposition. There 
are two particular interesting permutations: (/) the minimum 
permutation (0, 1, 2, . . . , ?7 — 1) whose canonical reduced de- 
composition is the identity, and (ii) the maximum permutation 
(n — 1, n — 2, . . . , 0) whose canonical reduced decomposition 
reads Si(s2Si)(s3S2Si) . . . {sn-i ■ ■ ■ si) , and it is the permu- 
tation with the largest number of simple transpositions in its 
canonical reduced decomposition. Then, a circuit that pro- 
duces the state corresponding to the maximum permutation 
needs as many gates as the number of simple transpositions 
of its canonical reduced decomposition, that is -2^^^^. The 
Laughlin state contains this maximum permutation, therefore, 
its quantum circuit must have, at least, "-^'^^^^ gates, that is 
precisely the number of gates in our proposal. 

It is also possible to analyze the way entanglement grows 
along the circuit. In order to do this, we calculate how much 
entanglement each gate of the circuit generates, that is, we 



determine the increase of the von Neumann entropy between 
the two parts of the system separated by a given gate. 
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The von Neumann entropy between k particles and the rest 
of the system is simply the sum of the contributions of those 
gates that are in the row which separates the system in k and 

n — k wires. These gates are vj.^ ' for k + I < n' < n and 
the entanglement entropy reads 
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This expression recovers in a clean way the same result as 
the one found in Ref. fisll . which was proven exact. Let us 
also remark that although each single particle is maximally 
entangled with the rest of the system, a subset of fc < n/2 
particles does not saturate the entropy. 



Sn,k < k log2 



(13) 



A experimental realization of our proposal will probably 
need to work on qubits rather than qudits. It is, then, necessary 
to find an efficient reduction of our algorithm to qubits. The 
easiest way to encode a qudit in terms of qubits is the binary 
basis. Then, an arbitrary single state \i) can be decomposed as 
\i) = \ir) ■ ■ ■ |«2)|ii) , where i = Y^l^i S'^^^ifc, r ^ log2 ?i 
is the number of bits needed to represent n, and ik ^ 0,1 
\/k = 1, . . . ,r. 

Now, we want to find the gates that acting on these 
qubits implement the W gates. Wij acts non trivially on 
the space spanned by the computational basis states = 
\ir, . . -,11, jr, ■ ■ ■ , ji) and \ = \jr, ■ ■ . , ji, ir, • ■ • and 
is the identity for the rest of states. Let us define W as the 
non trivial 2x2 sub matrix of W that acts on this subspace. 
According to Eq. W takes the form 



W 



Vp \/i 



(14) 



and it corresponds to the exponentiation W = exp (iOay) of 
the ay Pauli matrix for p — cos^ |. In order to implement an 
arbitrary Wij gate, we have to follow three steps: (/) first, 
we compare the binary expressions of . . . iij,. . . .ji and 
jr ■ ■ ■ jiir ■ ■ - H, and notice which bits are different. Then, 
we carry out a sequence of binary numbers, starting with ij 
and concluding with ji, such that adjacent members of the list 
differ in only one bit. These sequences are called Gray codes 
(ii) Next, we implement a quantum circuit performing 
a series of multi-qubit controlled gates that change the state 
ij according to the previous sequence. Each multi-qubit gate 
transforms the corresponding state of the sequence into the 
next one. These multi-qubits gates are carried out until it only 
remains a different bit between the last transformed state and 
ji. (Hi) At this point, we perform a controlled-M^ gate, or 
alternatively its complex conjugate W^, taking this different 
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FIG. 5: Implementation of the W35 gate in terms of multi-qubit con- 
trol gates. Each of these controlled operations can be decomposed in 
single qubit and CNOT gates. 

qubit as target. We will apply W or depending on the ini- 
tial state in which the single gate is performed in order to tune 
the behaviour of W according to its definition in Eq. (|6]l. Fi- 
nally, the reversed previous sequence of multi-qubit controlled 
gates is performed. 

This abstract construction can be illustrated with the ex- 
ample of the VF35 gate that acts non-trivially on the states 
|011 101) and jlOlOll). One possible sequence of Gray 
codes that connect Oil 101 and 101 Oil is 

Oil 101 

111 101 

10 1 10 1 (15) 

10 1 1 

10 1 Oil. 

From this, we can read its corresponding circuit, shown in 
Fig. |5] Notice that the first three multi-qubit controlled gates 
transforms |011 101) into |101 001). Next, the W gate is ap- 
plied to the fifth qubit, affecting only the states |101 001) and 
|101 Oil) due to the conditions on rest of qubits. Finally, we 
reverse the application of the multi-qubit controlled gates, en- 
suring that |101 001) gets swapped back with |011 101). It is 
important to point out that these multi-qubit controlled gates 
are not two-qubit gates, as, in principle, it would be suited. 
Nevertheless, it is known that these controlled operations can 
be performed by means of 0{r) single qubit and CNOT gate s 
iHlll which can be implemented experimentally lll8l ll9ll . 
Thus, if we consider that the number of M^-gates that we 



have to perform to implement our proposal for n qudits is 
n(n-i){27i-i) ^ number of single qubit and CNOT gates 
required by our circuit scales as O (n'^(log2 n)"^). 

There's a second way of encoding a qudit in the state 
that is, take n qubits, set the i-th one to |1), and then the rest 
to |0), i. e. \i) = |0 . . . 010 . . . 0). This unary encoding is 
less efficient than the previous one, since it requires n qubits 
compared to the logj n required before. Nevertheless, it has 
the advantage that allows us to implement the VF-gates using 
less gates. We just follow the previous three steps with the 
four qubits that have to be affected by the gate, that is, the 
i-th, j-th, {n + i)-th and [n + j)-th. We realize that only 3 
different multi controlled gates are required. In this case, then, 
the number of single qubit and CNOT gates required by our 
circuit would be O (n^) . 

Moreover, let us note that we can transform our antisym- 
metrization circuit into a symmetrization one by just chang- 
ing the signs of the definition of the W gates in Eq. (|6]l, 

i-e- ^/p\ii) T \/l --PIj*) ^ VpKj) ± VI - pli*)- 
Another possibility of performing the same symmetrization 
would be to invert the order of the input state of the circuit 
(jn — 1, . . . , 0) instead of |0, . . . , n — 1)) and to apply the 
gates Wn-i-i,n-i-j, instead of Wij, along the circuit. In 
both cases, the W gates always act in the positive combination 
and the final state obtained is fully symmetric in all possible 
permutations. 

Let us conclude with some comments on the generaliza- 
tion of our proposal to other values of m. The first obser- 
vation is that, if m > 1, the number of states that appear 
in the superposition of the Laughlin wave function is much 
larger than the simple permutations of the input single states. 
The corresponding quantum circuit, therefore, cannot be only 
composed of W^-gates. This will increase a lot the degrees 
of freedom of our elementary gates and, thus, its complexity. 
Though specific examples for low values of n and m can be 
found, a general scheme is still missing. 
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